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2 Gaussian Hopfield 

1. Introduction: Ising spins with a rotation symmetry 



In this paper we will illustrate the notions of chaotic size dependence, metastates and 
their dispersal, and the chaotic pairs of states scenario, introduced as a possible description 
of the low temperature spin glass phase [N,NS2,NS3,NS4,NS6], on a simple model which is 
similar to the two-state Hopfield model. The fact that the model has site disorder makes it 
more tractable than the commonly considered bond-disorder spin glass models. The main 
difference with the standard Hopfield model of neural networks, is that instead of two i.i.d. 
Bernoulli random variables the disorder is described by two i.i.d. Gaussian random variables 
at every site. As a consequence, in the thermodynamic limit we obtain the existence, for a 
"two-pattern" model, of uncountably many (instead of two times two) pure states for this 
model, due to the existence of a continuous (rotation) symmetry of the distribution of the 
random variables describing the disorder. In any finite volume, however, this symmetry is 
necessarily randomly broken in a given realization. Intuitively, this means that there are 
only two pure ground states, and the low temperature Gibbs state is close to the symmetric 
mixture of two, out of a possible continuum, of pure Gibbs states, due to the fluctuations in 
the disorder. 

The concepts we want to illustrate have their origin in the theory of spin- glasses. However, 
the most often considered spin-glass models, which have bond disorder, both in finite dimen- 
sion (the Edwards- Anderson models) and the equivalent neighbour (Sherrington-Kirkpatrick) 
model, have turned out to be so complicated to analyze, that up till now it has not been 
possible to check which of the possible scenarios for the spin-glass phase applies to them. 

We remind the reader that in the debate within the physics literature on the extreme sides 
there are the proposals of Fisher and Huse, [FH1,FH2,FH3,FH4] predicting the existence of 
only two pure states in any dimension higher or equal than 3, versus the proposal of Parisi 
and coworkers, in which an infinity of pure states is predicted [MPV, MPR]. This scenario has 
been claimed to apply down to the 3-dimensional Edwards-Anderson model. Intermediate 
scenarios have been discussed by [BF,NSl,NS2,NS3,NS4,NS5,NS6,N,vE]. 

Although of course lattice models with two pure states are common, our experience with 
models having an infinite number of pure states is a lot more limited. Therefore we hope that 
our discussion will be useful in illustrating various concepts, mostly introduced and studied 
in a systematic way by Newman and Stein (see in particular [N,NS2,NS3,NS4,NS6]), which 
have been introduced either in an abstract setting or via (in)formal arguments, by applying 
them to a concrete model. 
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The main idea in the approach of Newman and Stein is to classify the possible scenarios 
on the basis of first principles, using only general ergodic properties using the concept of 
"metastates" , i.e. probability distributions on the space of Gibbs measures (first introduced 
apparently in [AW]; see [N,NS2,NS3,Kul,Ku2,BG3] for more details, as well as applications 
of these concepts and extensions to equivalent neighbour or mean-field type models-to which 
our model also belongs). 

In this context, in one of their most recent papers [NS6], they conjectured that in a 
disordered lattice system, in any approximate decomposition of a finite volume Gibbs states 
into "pure states", the weights in this decomposition should be mostly concentrated on a 
single subset of states that are related by an exact symmetry of the system, while other states 
would appear with a weight that tends to zero as the volume tends to infinity. The particular 
subset chosen could of course be random and could depend strongly on the volume. Applied 
to the Ising spin glass situation, this argument would predict the chaotic pairs picture. 

Although a similar situation has been shown to occur in the usual Hopfield model with 
M = aN patterns if a is small in [BG3], we found it worthwhile to construct a simple model 
showing these features in order to see what is involved. 

Let us state the definitions of our variant of the Hopfield model and the main quantities 
of interest. Let <Sjv = {—1, +1}^ denote the set of functions a : {1, . . . , N} — > {—1, +1}, and 
the set S = {— 1, +1} N . We call a a spin configuration and denote by cr, the value of a at i. 
Let (f2, T, P) be an abstract probability space and let £f [u], i € N, (i = 1, 2, denote a family 
of i.i.d. standard Gaussian variables. We will write ^[uj] for the iV-dimensional vector whose 
ith component is given by [u]; such a vector is called a pattern. On the other hand, we will 
write £i[uj] for the two dimensional vector with the same components. When we write 
without indices, we consider it as a 2 x N matrix (its transpose will be denoted by £*). 

Throughout the paper, (•,•) denotes the scalar product, without indication of the space 
where its arguments lie. 

We define random maps m^[uj\(a) : Sn — ► [— 1,+1] (conventionally called overlap param- 
eters) through 




(1.1) 



i=l 
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The Hamiltonian is now defined as 

V 

m=i.2 (1.2) 

= -- 2-||"»jvH(<t)|| 2 , 

where || • || 2 denotes the ^2-norm in M 2 . 

Note that if we rewrite f'J = £ e ] = # cos(fl) + £ 2 sin(0) and £' 2 = £' ? = # sin(#) -£ 2 cos(#) 
the Hamiltonian has the same form in the primed variables. However, this transformation is 
a statistical symmetry, mapping one disorder realization of the model to another one, drawn 
from the same distribution, as opposed to for example the spin-flip symmetry which is an 
exact symmetry for any given realization of the disorder. 

Through this Hamiltonian, finite volume Gibbs measures on Sn are defined by 

p -j3H N [Lj]{a) 

and the induced distribution of the overlap parameters 

QnA u \ = A*iV,^M ° rriNluj}' 1 . (1.4) 

The normalizing factor in (1.3), called the partition function, is explicitly given by 

Z NtP [u] = 2~ N e-f HN ^ = E a f^HW. ( L5 ) 

We are mainly interested in the concentration behaviour of Qn,p as N — > oo. It will be con- 
venient to do this by considering the auxiliary measure Qn,/3 = Qn,/3 *A/2(0, ^1) obtained 
by a convolution with a Gaussian measure, its so-called Hubbard-Stratonovich transform. 
Since, for N large, A/^O, ^1) converges rapidly to the Dirac measure at zero, the two mea- 
sures have asymptotically the same properties. For details see e.g. [BGP]. Qn,p is absolutely 
continuous with respect to Lebesgue measure on M 2 and has the density 

-PN^ N ^[w\{z) 

7 r„i . (1-6) 

where 3>jv,/3 is given by 

1 1 N 

$jv^M(ct) = -\\z\\ 2 2 - — J^ln cosh/3feH,z). (1.7) 
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As usual in mean-field models, we construct the extremal Gibbs measures by tilting the 
Hamiltonian (1.2) with an external magnetic field (for a general discussion on the issue of 
limiting Gibbs states in mean field models, see [BG1], Sect. 2.4 or [BG3], Sect. 2). That is, 
we define a more general Hamiltonian 

H%[u]{a) = -^\\ mN [u]{a)\\l - N(h,m), (1.8) 

where h = (6cos(i9), 6sin(i9)) £ M 2 . The corresponding measures on the spins and on R 2 are 
denoted by n% p[uj] and Q% p[u>], respectively. We then take the limits lim&^o limjv->oo) for 
all values of ■& € [0, 2tt). We distinguish the measures constructed from this Hamiltonian by 
an additional superscript h. 

We are now able to give a precise formulation of our main results. 

Theorem 1: Let h = (6cos#, bsmfi). Then 

lim lim Q% g = 6r r * C ostf,r* sin*)) (1-9) 

where r* is the largest solution of the equation 

dxe~^xtanh(/3xr*). (1.10) 



J 



Theorem 1 shows that there is an uncountable number of extremal limiting induced mea- 
sures, indexed by the circle. The following Corollary shows that to each of them corresponds 
a distinct limiting Gibbs measure on the spins. 

Corollary 2: For any finite set IcN, and F-almost all u, 

*WM(K = sj}) ee UmJim^>](K = = II 2c0Bh(j9( g iM>m)r (1-H) 

w/iere m = (r* cos(#), r* sin(i?)), and r* as in (1.10). 

In Theorem 1 and Corollary 2 convergence is almost sure due to the presence of the tilting 
field. The situation changes if we set b = first and take the infinite volume limit later. 

Theorem 3: Let Qn,/3 as in (1.4) and m = m(i?) = (r* cos$,r* sin#) ; where d £ [0, tt) is 
a uniformly distributed random variable. Then 

v 1 1 

Qn,0 ~<*m(0) + ^-m{&) = QooA m \- ( L12 ) 
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Furthermore, the (induced) AW-metastate is the image of the uniform distribution offi under 
the measure-valued map # t— > Qoo,p[m{'&)\. 

Corollary 4: Let I C N be finite. Then the following holds: 

(i) Let {gi} i€ i be a family of i.i.d. random variables, distributed as jV(0, r*). Then 

hm = SI ) % \ g + \ g 2^h^"- ( L13 ) 

(ii) The AW-metastate is the image of the uniform distribution on d under the measure-valued 
map ■d i-> Aioo,mwM w/iere 



J e /3s»(£i[w],m) J g-^Si($»H,m) 

Moo,/3,mM = g II 2cosh/3feM,m) + 2 n 2cosh/3(6^,m)- (L1 ' ) 



Statement (ii) of Corollary 4 motivates the notion of metastates. Whereas on the level of 
the induced measures Qn,/3 one cannot see any influence by the conditioning, this is clearly 
the case on the level of the Gibbs measures on the spins. 

The remainder of this paper is mainly devoted to the proofs of the two theorems (the 
corollaries are standard consequences (see e.g. [BGP1] or [BG3] for proofs of analogous 
statements in more complicated situation) and will not be given) is organized as follows. In 
Section 2 we prove the necessary concentration estimates on the measures Qn,/3- This will 
yield immediately Theorem 1. In the case h = we will show that the measure concentrates 
near the absolute minima of some random process, and in Section 3 we will analyse the 
properties of these minima. In particular we will prove that these converge in distribution to 
one-point sets. This will allow us to prove Theorem 3. In Section 4 we discuss some further 
consequences on the chaotic volume dependence, the empirical metastate and the superstate. 

Remark: We consider the case of two patterns here in order to keep technicalities to a 
minimum. All our results can be extended without any novel difficulties to the case of any 
fixed finite number, M, of Gaussian patterns. In that case the set of extremal Gibbs measures 
will be indexed by the sphere in W M and the metastate will be supported on pairs of mirror 
images on this sphere, with the position being uniformly distributed. Thus nothing really 
new will happen. The situation when the number of patterns grows with the volume may be 
more interesting and work in this direction is in progress. 
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2. Concentration 

In this section we show the concentration properties of the measures Qn for large (3. These 
imply the same concentration results for the measures Qn by standard arguments that have 
been developed in much more complicated situations, see e.g. [BG2]. The estimates presented 
here are mostly similar, and often much simpler, to those that can be found e.g. in [BG2], 
but we decided to present some parts in detail where some care is required. 

We start with the more delicate case h = that will be relevant for the proof of Theorem 3 
(which will be given at the end of Section 3). We are interested in the concentration behaviour 
of the measures Qn,/3- The following two lemmata each give a partial answer. The first one 
asserts that Qn,/3 is concentrated exponentially about a circle around the origin, whereas the 
second one tells us that even on this circle, only a small part really contributes to the total 
mass. 

Lemma 2.1: Let {£f }ieN,/*=i,2 be i.i.d. standard Gaussian variables, and define <&N,p{z) 
as 

1 \ N 

<S>n{z) = -HI - — J> cosh/3fe,z). (2.1) 

Let furthermore Sjy = iV -1 / 10 . Then there exist strictly positive constants K,K', m,m' 
such that (r* is the largest solution in (1.10)) 

f -f3N® N (z) d 
J\ M\-r-\>S N < KN ~ 

f\ M -r-\<6„e- f " f * t,Wd *- 

on a set of ' P '-measure at least 1 — K'e~ K N . 

The second result needs an additional definition. Let 

1 N 

g N {$) = — V In cosh(/Jr*C< cos(tf - (2.3) 

where (Q, ipi) are the polar coordinates of the two dimensional vector £j. 

Lemma 2.2: Assume the hypotheses of Lemma 2.1. Let = iV -1 / 25 . Then there exist 
strictly positive constants K\ , K 2 , C\ , C2 such that on a set oj "P-measure at least 

»7~l/25 

1-K ie - N (2.4) 
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the following bound holds, 

f A , e-^^^dz 



f e -f3N* N (z) dz 
J An 



<C ie - N2/ \ (2.5) 



where 



A N = {(r,0) G 1R+ x [0, 2tt)| \r - r*\ < S N , g N (#) - ming N (#) < a N \, 

(2.6) 

A' N = |(r,0) G M.q x [0, 2-7r) 1 1 r — r* | < fo, 5iV (tf) - rnin^O?) > a N j. 

Combining these two lemmata and using the Borel-Cantelli lemma, we get immediately the 
following result. 

Proposition 2.3: Assume the hypotheses of Lemma 2.1. Then there exist strictly positive 
constants K,K',m, such that 



r e -f3N* N (z) dz 
J An 

where A^ is as in Lemma 2.2. 



>Ke~ K ' Nm , t.o.tnN 



0, (2.7) 



To see why the preceding results should be expected, we must consider the function &n,p- 
Note that the expectation of this function, 

E$ N (z) = -||z||| - — E In cosh/3(£i,*). (2.8) 

depends only on the modulus of its argument. It is useful to observe that if z = (r cos 6, r sin 0), 
we can represent E&n(z) as 

E® N (z) = K 2 -E^E c ln cosh(/3rCcos(^))d99 (2.9) 

where (, (ft are the representation of the polar decomposition of a two dimensional normal 
vector, i.e. £ is distributed with density xe~ x ' ' I 2 on R+ , and ip uniformly on the circle [0, 2ir). 

From this it follows that E<&jv(z) takes its minimum on the circle with radius r*(fi), where 
r* is defined in Theorem 1. It is easy to verify that there is < (3* < oo, such that r*{0) > 
if and only if (3 > f3* . 

It is also straightforward to check that E $ is sufficiently smooth to guarantee that it is 
bounded from above by a quadratic function (of ||z||) in some neighbourhood containing r* . 
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Proof of 2.1: We start with the numerator. We decompose the domain of integration into 
an "inner" part T,and an "outer" part O: 



{zeR 2 : \\\z\\ - r*\ > 5} = {z G R 2 : - r* > 5} 

u{z£R 2 : -r* < -5} = Ol)l. 

Consider the integral on O. We write it as 



(2.10) 



f e -N$ N (z) dz = I e -/3NE<t> N (z) e -/3N{<f> N (z)-E<Z> N (z)) ^.11) 
JO Jo 

and observe that E&n can be bounded below by a quadratic function C(||z|| — r*) 2 . We 
are left with the task of estimating the term &n{z) — K^n(z). This is accomplished by the 
following Lemma. 



Lemma 2.4: Let /jv(z) = Z)ili ln cosn an< ^ 

= {zeR 2 : \\z\\ > r* + 5} . 



(2.12) 



Then, for 5 small enough, such that 5 2 /16 < 5/2y/2, there exist strictly positive constants 
C\,C-2,K\, K 2 such that 



C 



*\2 



snp\f N (z)-Ef N (z)\>-(\\z\\-r*) 
zeo 2 



< Kl e~ K ' N + C^e-^^N- 1 *. (2.13) 



Proof: Define /n(z) = /n(z) — E/jv(z). The left-hand side of (2.13) is bounded from above 
by 



< 



C 



sup |/*(;0|>- r (||s , ||-r* 



z'€W r no 



+ 



sup sup \/n{z) - In(z')\ > 

z'e\VrnAz€Br.(z') 



c 



\z'\\-r* 



(2.14) 



where W r is the grid with spacing r in I 2 , and z' £ W r is chosen such that < \\z\\ — \\z'\\ < 
V2r. 

The argument of the second term can be uniformly bounded. Using e.g. Lemma 6.10 of 
[BG1], we get that 

IM^)-M^)l<lk-^ll2PH 1/2 , (2.15) 

where A is the matrix (l/iV)£ T £. Similarly, 

\Ef N (z)-Ef N (z')\ < Wz-z'UnAWf/ 2 . (2.16) 
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Now, a trivial computation shows that 

E||A|| < 1 + C/y/N 



(2.17) 



and using (for instance) the same argument as in Section 4 of [BG1], but replacing Talagrand's 
concentration estimate for bounded r.v.'s by the standard Gaussian concentration inequality 
(see e.g. [LT], Ch. 1), one shows easily that 



\A\\ - II > x\ < Ce~ Nx2 / c . 



(2.18) 



Therefore, 



C 



sup sup \f N ( z )-f N ( z ')\>^\\ z '\\-r*^ 



z'ew r nA z eB r (z>) 



< 



< 



C 



r(||A|| 1 /2 + (E p||)i/2 ) >_ ( ||^||_ 



*\2 



(2.19) 



(AH 1 / 2 + 2) > 



CS_ 

4r 



2 1 



Choosing the grid parameter r such that r < C5 2 /16 the right-hand side of (2.19) is bounded 
by P[||A|| > 4] < Ce~ 9N/c This takes care of the second term in (2.14). Let us now treat 
the first term. The probability that the supremum over all lattice points of some function 
exceeds some given value is transformed into a summable series of probabilities that at each 
lattice point the function is greater than this value. More precisely, we have 



sup \f N (z')\>j 
U'evvvno 4 



\z'\\ -r*f 



< £ P \f N ( z >)\>^ { \\ z >\\- r *f 
z'ew r no L 

< ^ e -KC 2 (\\z'\\-r') 4 N^ 



(2.20) 



z'ew r no 



by Chebyshev's inequality. Then 

-KC 2 (\\z'\\-r') 4 N = r -2 Y 



z'ew r nO 



r 2 e -KC 2 (\\z'\\-r") N 



z'€W r no 



< 



r->f 
Jr.- 



e -KC 2 (\\z>\\-r*) N dz 



< r - 2 e- K ^ 5iN 



B {r"+8-V2r) 

[ e- K ^'W-r*Y» dz 

JVL 2 \B {r»+&/2) 

ze~ KzA dz 

6/2 



(2.21) 



4 / 

Js 

< K'r- z e-*^°' n N't, 
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where K 1 stands for an upper bound for the integral, which is independent of N (assuming 
5 > 2y/2r). Combining this and (2.19), and choosing S small enough such that CS 2 /16 < 
6/(2y/2) concludes the proof of Lemma 2.4. 

Therefore, on a set of measure at least 1 — C\e~ c,lN ^ , the integral (2.11) can be bounded 

by 

f e -pNV*,,{z) e -pN{* N {z)-E* N {z)) fa < f e -f3N<f(\\z\\-r*) 2 dz 

Jo Jo 

/oo 
re -p NC{ r-r* f dz 

t+S oo ( 2 - 22 ) 



sup|/jv(2) -Ef N (z)\ > e 
zeA 



= 2vr^-e-^ 2 . 

pNC 

We now turn to the integral on the "inner" part I. Again, we have to control the term 

$ N (z)-E* N (z). (2.23) 
Since X is compact, we can do this uniformly by using the following lemma. 

Lemma 2.5: Let f N (z) = l/((3N) In cosh/?(&,z) and A C M 2 a bounded set. Then 
there exist strictly positive constants K\ , K 2 , C\ , C2 such that 

The proof is similar (if not simpler) to the proof of Lemma 2.4 and is left to the reader. 
Lemma 2.5 implies that 

f e -f3N^ N(z ) dz < e sN e -f3NEHr*) f e -0NE* N {z) ^ 

Jx Jx (2.25) 

using the fact that E$jv(||,z||) — E<£(r*) can be bounded uniformly on 1 by its value for 
\\z\\ = r* — 5. 

Finally, the denominator in (2.2) can be bounded from below, using the second order 
Taylor expansion with remainder of E $jv(||z||) 

f e-^^ z Uz 

J\ \\z\\-r*\<5 



> e -/3VE*(r*) f e - NC (\\z\\-r*) 2 -NC'{\\z\\-r*) 3 -Ne dz ^.26) 

J\\\z\\-V\<6 
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on a set of measure at least 1 — Ke~ KN — Ce~ 2 e~ CNs (this error term can be estimated by 
Lemma 2.5). Collecting (2.22), (2.25) and (2.26), we get that on a set of measure exponentially 
close to one, 



r -0N3> N { Z ) j 

J\\\z\\-r*\>5 e 

J | ||z|| — r* <0 



< Me £ ' 3Jv e^ c '' 53 (27r)- 1 /3iVC (l - te"^ 2 ) 



(3NC \ (2.27) 

MKe -flN(C6>-2e-C>6*)p N ^ _ ^-A^" 1 

+ Mif'e-^(^ S - £ - c ' i3 )JV (l - ,5e-' 3JVC ' 52 )" 1 . 



Now let us choose 5n = TV io , £jv = N * ; then (2.27) gives 

r -pN$ N (z) dz 4 j 1 

Jj H^ll — r-- |><5jy < jyj-g-jy e -l3N5 (C-2N~?<5 -C' N'TS) 

f lM _ r , l<SN e-W*»(*)dz ~ (2.28) 

+ M^7Ve- iV * ( 4- JV ~*- c ' iV ~*), 

on a set which is exponentially close (in N) to 1. This concludes the proof of Lemma 2.1. 

We now turn to the proof of Lemma 2.2 which is a little more delicate than the previous 
one. 

Proof of 2.2: Let us write 1(B) for the integral J B e~ /3N ^ N( - z ^ dz. We will prove the 
concentration behaviour by a strategy similar to the one used in Lemma 2.1. Namely we 
replace the function $at by its expectation E$at and control the error. 

Write the fluctuation term <E> jv — E at as 



1 N 

$ N (z) -E<f> N (z) = — ^{ln cosh/3(6, z)-E\n cosh/?(&, z)} 
1 * 

= ]^Y^ ln cosh P(& z ) - ln cosh ^ ^ 

- E In cosh /?(&, z) + E ln cosh /?(&, z')} 
1 N 

Now choose z' such that z' = z'(z) = Xz, A > 0, and ||z'|| = r* (i.e. z' is the projection of z 
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onto S' 1 (r*)). Define the two functions 



1 N 

h N{z) = -7^y2i ln cosh P(tii,z) - In cosh/3(&,z') 
v N ,=1 

- E In cosh/?(&,z) +E In cosh /?(&, z')}, 



(2.30) 



with z' defined as above, and 

N 



g N (z) = -J= V{ln cosh /?(&,*) -E In cosh/?(&, z)}. (2.31) 



i=l 



Then the fluctuation term takes the form 

N(* N (z) -E$ N (z)) = ^-{h N (z) - g N {z')). (2.32) 

It is the term that determines the concentration behaviour of the measure. To see this 
we first bound the term uniformly on the "annulus of concentration" U A' N . We have 
the following result. 

Lemma 2.6: Let {£i}ieN be i.i.d. Gaussian variables with mean zero and variance one. 
Let tiN be as in (2.30), and A' N as in (2.6). Then for any e > 0, 



sup \hjs[(z)\ > £ 



zeA N L)A' N 



<KN 2 e- Nl/1 °^- KN ' 1/10 \ (2.33) 



Proof: Let us write 

fi(z) = lii cosh /?(&,*), (2.34) 

and 

fi = In cosh/?(&,z) -E In cosh/3(&, z). (2.35) 

We also keep the notation z 1 = z'(z) defined above. Introduce a polar grid WV in M 2 , i.e. a 
discrete set of points Xij whose polar coordinates are given by (pi,ctj) £ M + x [0, 2n), such 
that Ajya = |q;, — ctj\ = KN~ 1 / 2 and Afjp = \pi — pj\ = KN -1 / 2 , for some appropriate 
constant K. Note that for any point z in a bounded domain i C R 2 , the distance to the 
closest grid point is less than K 1 N~ x l 2 . 

For any z 6 M 2 , define x = x(z) € Wat to be the grid point closest to z, and y = y(z) G Wat 
the grid point closest to z' = z'(z). One can easily convince oneself, that x' = y' , i.e. the 
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two points x and y lie on the same ray starting at the origin. Then we can decompose the 
function }in{z) as 



h N (z) 



1 N 



i=l 
N 



^ (z) - fi{x)} + -L - My)} 

1 v i=l 



(2.36) 



i=l 



1 " 

+ 7^E^(y)-/^')}- 

V i=l 



Denote by h(z, x), h(x, y), h{y, z') respectively the first, second and third sum on the right- 
hand side of (2.36). We can then write (let An = Am UA' n , the "annulus of concentration") 



sup |/ijv(z)| > £ 

z€An 



sup \h(z,x) +I 2 (x,y) +I 3 (y,z')\ > e 



< 



sup 



x£WnDAn zeB KN _ 1/2 (x) 



sup \Ii(z, x)\ > 



+ . 



+ : 



sup sup |/2(x,y)| > | 
xew N nA N v^ w N n - A N <J 

y' =x' 



(2.37) 



sup 



yew N nA N z'eB KN _ 1/2 { y ) 



sup \I 3 (y,z')\ > 



The first and the third term (they are equal) can be uniformly bounded by an estimate 
analogous to the proof of Lemma 2.2. In fact, for any u, v, we have 



1 N 

^= £{£(«) -/*(<;)} 



N 1 



< v / ]v/3(Pir /2 + (Ep||) 1 / 2 )||u- 



V 2- 



(2.38) 



Now, if || u — v || 2 < Ae'N 1 / 2 //3, we have the following exponential bound. 



N 



— J2{fi(u) - Uv)}\ > s' 

i=l 



< 



[A^/ 2 + (E\\A\\) 1 / 2 > 



e'N- 1 / 2 
3\\u - v\\ 2 



(2.39) 



< P[P|| > 4] < Ke 



-KN 



Thus we get for the first term in (2.37), 



sup sup \h(z,x)\ > - 

x€W N nA N zeB KN _ 1/2 (x) ^ 



xeW N nAN 

< ^ P[||A|| > 4] < i^iV 1 / 10 ^-^ 



sup \h(z,x)\ > - 



(2.40) 



-KN 



x€Wn<~)An 
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since we know that \\x — z\\ = K' 'iV -1 / 2 , by the remark preceding (2.36), and the number of 
grid points in An is bounded by N x 5^ 1 times some constant. The same estimate is valid for 
the term containing I 3 (since they are equal). 

Let us now consider the term containing I 2 . We know that \\x — y\\ < 25 n, since those two 
points are supposed to lie on the same "ray" . Again, we can turn the supremum into a sum, 



sup sup \h{x,y)\ > - 
x€WnI~\An v£W n da n o 
«'=*' 



<£p[|/ 3 (y,z')l>f], (2.41) 

x,y 



where x,y on the right-hand side satisfy the same conditions as on the left-hand side. By 
Chebyshev's inequality, we get that for any u, v 

N 

< inf e- S£,VW E 



.i=i 



s>0 



inf e~ S£ 

s>0 



N 



(2.42) 



JjEe»{/*(«)-/i(")}. 
i=i 



Now we use the series expansion of the exponential function, the fact that the exponent in 
the right-hand side of (2.42) is a centered random variable, and some obvious inequalities for 
each term of the expansion, to get 

[(/,(«) - fiWfe'MM-ftW] } . (2.43) 



(2.44) 



]g e s{/»H-/iW} < J i _|_ 



To evaluate the expectation term, we use the inequality 

\fi(u)-fi(v)\ <P\(Zi,u-v)\. 
Then the expectation term in (2.42) is bounded by 



E 



(fi(u)-fi(v)) 2 e a MM-f*W 



< (E [(/ 4 («)-/i(t;)) 4 ]) 4 (Ee 2s ^-f^) 

< 4 (E [(fi(u) - h{v)f]Y (Ee 2s ^~ f ^ 

,sE\fi(u)-fi(,v)\ 



(2.45) 



x e° 



where the first inequality follows by Cauchy-Schwarz, and the second one is a consequence of 
the inequality (a + b) 2 < 2(a 2 + b 2 ) (applied twice to the first factor), respectively the trivial 
fact that |a — b\ < \a\ + \b\. All quantities in (2.45) can be bounded easily using (2.44). One 
gets (by calculating explicit Gaussian integrals) 



E [(fi(u) - fi(v)) 4 ] = 3\\u - 



12! 



(2.46) 
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]g e 2s|/i(u)-/i(t))| < 2 e 2s2 W u - v \\2 ^ 
e s¥,\fi(u)-fi(v)\ < e sy/2/n\\u-v\\ 2 _ 

Inserting (2.46)-(2.48) into (2.45), gives 
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(2.47) 
(2.48) 



-E 



(fi(u) - fi(v)) 2 e slMu) - f ' {v)l < 2VQs 2 \\u-v\\ 2 2 e 2s ll«-«ll 9 W 2 /*ll«-«lb. ( 2 .49) 



We use the above bound (2.49) in (2.42), together with the inequality 1 + x < e x , and the 
fact that \\x — y\\2 < <5jv = KN -1 / 10 . We thus get the following estimate 

N 1 

< inf e -^'VN+K a a N^'e a ' aN ~ 1/B +^ N ~ 1/10 



1=1 



s>0 



(2.50) 



Choosing s = N 2 / 5 , this gives 



JV 



Y,{Mu)-Mv)}>VNe' 



i=l 



< Ke~ Nl/10( - e '~ KN ~ 1/W) 



(2.51) 



The same bound applies to 



N 



i=l 



(2.52) 



Inserting (2.51) and (2.52) into the left-hand side of (2.41) gives 



sup sup \I 2 (x,y)\ > e' 

xew N nA N y^N n - A N 



< KN l/2 N l/10 e -N^o {£ ,_ K , N -mo^ (2 ^ 3) 



since the number of terms in the sum does not exceed a constant times N 1 / 2 (the number of 
allowed x) times iV 1 / 10 (the number of allowed y). Using (2.40) and (2.53), (2.37) gives 



sup \}In(z)\ > £ 



This concludes the proof of Lemma 2.6. C> 



< KN 2 e- K ' Nl/1 ° s . 



(2.54) 



Note that we can choose e as a function of N, and still get an exponential bound. For 
example, choose e = = (In N) 2 N~ l / 2Q . Lemma 2.6 then reads 

Lemma 2.7: Let {£i}ien be i.i.d. Gaussian variables with mean zero and variance one. 
Let Hn be as in (2.30), and An , A' N as in (2.6). Then, 



sup 

IM*)| > iV" 1/20 (lniV) 2 

zeA N UA' N 



< KN 2 e -N^ { ^Nf-K'N-^)_ (2 _ 55) 
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Furthermore, 



sup 

IM*)| > iV" 1/2 °(lniV) 2 , i.o. in N 

zeA N UA' N 



= 0. 



(2.56) 



Proof: The first statement (equation (2.55)) is a straightforward consequence of Lemma 
2.6. Equation (2.56) then follows by the first Borel-Cantelli Lemma. 

Let us now estimate the integral I(A'n). We get explicitly, using the bound on from 
Lemma 2.6, 



e -PN* N {z) dz 



e -(3NE<i> N (z) e -^Nh N (z) e -y/Ng N (z'(z)) ^ 



</ 

J\r-r*\<5r, 

L 



re -0NE<S> N (r*) e y/Ne 



'Siv(i?) — min g N >a N 
= 2e -f3NE^ N (n e VNe I rdr 

\r — r* \<8n 



I 

J OA 



(2.57) 



x / 

'SiV-(^)- min 9N>a N 
< 4 e -/3NE<S> N (r*) e VNe r *fi 



Thus, 



AT 

x 2ne~^ r ^ aN e~^~^ min9N 

e -f3NZ N (z) dz < Xe -/3JVE$ w (r*) e ViV £(5Arr * e -VJVa w _ 



(2.58) 



We now turn to the integral I(An)- Using standard estimates for Gaussian integrals, a 
quadratic upper bound of about its minima, and the fact that E$(||z||) can be bounded 
from above by a quadratic function in some neighbourhood containing r* , we get 

f e -PN^ N (z) dz > e -0NE<S> N (r*) e -VNe f r( ,-PNC' (r-r*f dr 

J A N J\r-r*\<5 N 



J gAr(i?)-min g N <a N 



v Arc / 



1/2 



(2.59) 



(1-e 



7T 



1/2 



(1-e 



-\/~N K'aN 



)■ 



We get finally for the ratio I(A' N )/I(A N ) 



I(A N )- K r*-8 N iy 



(2.60) 
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Lemma 2.7 allows us to choose e = e{N) = iV _1 / 20 (m TV) 2 . Inserting this choice, together 
with a N = iV" 1 / 25 , into (2.60), gives 

I {An) 

This statement is true for all w £ fi, for which Lemma 2.6 respectively 2.7 holds, that is on 
a set of P-measure at least KN 2 e~ Nl/2 °^ ln Nf-K'N- 1 / 20 ) _ Thig proves Lemma 2.2. 

Let us now turn to the proof of Theorem 1. We again state first a result about the 
concentration of the induced measure Q%,f3- 

Proposition 2.8: Let {£f }ieN lAt =i,2 be i.i.d. standard Gaussian variables, and define 

1 1 N 

*nA z ) = 2 l|z| ' 2 " Jn cosh ^ z + h )- ( 2 - 62 ) 

Let furthermore Sn = -/V -1 / 5 . Then there exist strictly positive constants K, K', m such that 

^ J||z-f h ||<<5jv C ' UZ J 

where r h is the unique minimum of the function 

®*nA z ) = 1^ ~ l Eln cosh/3(6,^ + h). (2.64) 



Proof: Let us decompose <& N a in the usual way 

^{z) = E<&* j/3 (z) + * h N j{z) - m h N #{z). (2.65) 

We first treat the denominator appearing in (2.63). E<3?^ ^ can be bounded from below by 
some quadratic function C\\z — r h \\2 on the set \\z — r h \\ >5n > 0. The fluctuation term can 
be controlled by the following analogue of Lemma 2.4. 

Lemma 2.9: Let /n = Yl^=i m cosh /3(^i , z + h). Then for 5 small enough, sucht that 
C5 2 /80 < 5/2, there exist strictly positive constants C\ , Ci , K\ , such that 



Pn 



C 

-f h \\ 2 >5 2 



sup \f N {z)-Ef N {z)\ > ^r\\z-r h \\l 



(2.66) 



< Kie -K 2 N + CiN l/2 § -2 e -C 2 



N 
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Proof: The proof is completely analogous to the proof of Lemma 2.4, and is left to the 
reader. 



Therefore, with probability greater than 1 — pjy, sup($>% a — EQ^^Jz) does not exceed 
one half of the lower bound of the deterministic part, which implies that 

e -f)N** ift (z) dz < e -/9JVE** >/9 (f*) f e -p N c llz _ fhlll ^ 

\\z-f h \\>S N J\\z-f h \\>5 N (2.67) 

< e -f3NE$ h N ^(? h ) e -f3Nq. S 2 N K^ 



L 



We now turn to the denominator in (2.63). The probability that the fluctuation term 
exceeds an e > is bounded by Lemma 2.5: 

q N = ^\ sup \f N (z)-Ef N (z)\>e\<K ie - K * N + C\e- 2 e- C262N . (2.68) 

L \\z-f h \\<5 N 

Using the Taylor expansion of ^(z) about f h up to order 2, with an error term of order 
3, we get that with probability higher than 1 — q^, 

[ e -l3N^ ]P (z) dz y e -f3N(^ h NJ3 (f h +C''S 3 N +e))f g - /3NC' \\z-f h f 2 dz 

J\\z-f h \\<6 N J\\z-f h \\<6 N (2.69) 

> e -pN(*%,p(r h +C"S a N +e)) KN -l/2( 1 _ e ~f3N^-8 2 N ^ 



Combining (2.67) and (2.69) gives 



J\\z-rh\\>S N e UZ / ^„-0N(£-5l, 



J ||z— r || <ojv 

with probability greater than 1 — (g^ +pn)- Choosing 5n = -/V -1 / 5 , e = iV -1 / 5 , implies 
that Y2n(pn + Qn) < °°- Applying the Borel-Cantelli Lemma then gives the statement of 
Proposition 2.8. 

Theorem 1 is now obvious: 

Proof of Theorem 1: Let / be a bounded continuous function. Then 

QnAV = f(.r h Q h N A*{\\*-i*\\<6»}) + Q h NAV(r h ~ /)I{||*-f*||<^}) 

+ QNAf^{\\z-r h \\>SN})- 

Taking the limit N | oo, we can replace Q% p by Q% p and use Proposition 2.8. Since / is 
bounded, the third term on the right-hand side of (2.71) converges to zero, and since it is 
continuous, the second term also vanishes too. These statements are true P-a.s. Finally we 
let b = \\h\\ 2 — > 0. Again by continuity of /, f(f h ) — ► /(r*(cosi?,sin#)). This proves the 
Theorem. 00 



(2.71) 
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3. Uniqueness of extrema of certain gaussian processes. 

In the previous chapter we have seen that the measures Qn,/3 concentrate on a circle of 
radius r* at the places where the random function 5jv(i?) takes its minimum. In this section 
we will show that these sets degenerate to a single point, a.s. in the limit N f oo. To do so 
we first prove a uniqueness theorem for the absolute minimum of a certain class of strongly- 
correlated Gaussian processes. Then we show convergence in distribution of <7jv($) to such a 
process and finally we show that this implies also the desired convergence in distribution of 
our measures. We begin with the following general result. 

Proposition 3.1: Suppose is a real stationary Gaussian process which is periodic 
with period T. Suppose furthermore that its covariance function r(s,t) = r(s — t) is even, 6 
C°°[0, T], and r(r) is less than r(0) for all r 6 (0, T). Then there exists an equivalent process 
rj(t) having almost surely infinitely differentiable sample paths. Moreover, the probability that 
there exist two or more maxima with equal height in [0, T) is zero. 

Proof: Without restricting the generality, we can assume that E [x(t)] = and re = 
E[x(*) 2 ] = l- 

By its continuity properties, r(r) can be expanded about the origin as 

r(r) = l-^ + 0(r% (3.1) 

The first assertion then follows from the following result due to Cramer and Leadbetter (see 
[CL]), chapter 9.2). 

Lemma 3.2: Suppose that for some a > 3, 

r(T) = 1 - blS + Q ,3.2) 

where A2 is a constant. Then there exists a process rj(t) equivalent to and possessing, 
with probability one, a continuous derivative rf(t). 

Proof: See Cramer /Leadbetter [CL]. 

It is easily checked that by (3.1), r(r) satisfies the condition (3.2) in Theorem 3.2, which 
proves the statements about continuity and existence of a continuous derivative. 

Consider now the process x'(£). Its covariance function f(r) is given by f(r) = —r"(r) 
(see for example Leadbetter et al. [LLR], p. 161, chapter 7.6). Then it can be expanded 
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about the origin as 



r(r)=A 2 -^r 2 + 0(r 4 ). 



(3.3) 



Then f(r) also verifies condition (3.2) in Theorem 3.2. Repeating this argument implies, 
together with the Borel-Cantelli Lemma, that there exists an equivalent process rj(t) having, 
with probability one, infinitely differentiable sample paths. 

From now on, we assume that X (t) itself has the above continuity properties. We want to 
find the probability that there are not two maxima with equal height in [0,T), i.e. 



P [3s, teTxT:\s-t\^kT, \ X (t) - X (s)\ = 0, | X '(£)| = \ X '(s)\ = 0] = 0. 
We first show that for any $ > 0, 



3s, t € T x T : 



kT 



>^lx(0-x(s)l = o,|x'(t)| = |x , (5)| = o 







(3.4) 



(3.5) 



Let us choose a collection of grid points £j € T, separated by some distance e > 0. By the 
continuity properties, X and X ' are Lipschitz-continuous with a.s. -finite constants Co, C\. 
Consider the set ricCri such that Co and C\ are bounded by some number C > 0. Then, by 
Lipschitz-continuity, x'W = 0, £ € [£j,£j+i) implies that (for some x G [ti,t]) 



\X'(U)\ < Ce. 



Similarly, \ X (t) ~ x( s )\ = implies 



\ X (t l )-x(t J )\<2Ce 



(3.6) 



(3.7) 



where t — ti < e, s — tj < e. Then we can estimate the probability of the event in (3.5) (on 
f2) by 



3s, t £ T x T : 



kT 



>$,\x(t)-x(s)\ = 0,\x'(t)\ = \ X '(s)\=0 



< 



3U,t 3 : \kT-\s-t\ \ >^| x ( ti )- x ( t .)| <2Ce,\ X '(t t )\ <Ce, 



(3.8) 



lx'(^)l < Ce 



Let us denote the event appearing on the left-hand side of (3.8) by A$, and the event appear- 
ing on the right-hand side by B# t£ . The probability P[^ ;£ ] can be estimated by the standard 
bound 



P [B*, e ] < Yl F MU) - X(h)\ < 2Ce, \ X '(U)\ < Ce, < Ce] . (3.9) 

\kT-\U-t j \\># 
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Now, for any fixed 

(x(*0-x(ti),x'(tO,x , (^)) (3-10) 

is a Gaussian vector, and due to the condition on |ij — and the assumption concerning 
r(r), its distribution is non-degenerate. Therefore, each term in the sum on the right-hand 
side of (3.9) can be bounded by 



' [\X(U) ~ X(h)\ < 2Ce, \ X '(U)\ < Ce, \ X '(t 3 )\ < Ce] < Ke 3 C 3 (27ra 



1,3 j 



(3-11) 



where cr^j is the determinant of the non-degenerate covariance matrix of the random vector 
(3.10). Since the ti,tj are chosen in a compact set, this quantity can be bounded uniformly 
in i, j. We thus get 



' l\x(U) - X(tj)\ < 2Ce, \ X '(U)\ < Ce, \x'{t 3 )\ < Ce] < K(d)e 3 C 3 . 



(3.12) 



Finally, the number of allowed pairs in the sum in equation (3.9) does not exceed T 2 e 2 , 
which implies that 



n 



c 



<K{$)T z e- z e 6 + 



(3.13) 



keeping track of the set £l c c on which the above estimates are not valid. Now choose C 
C(e) = o(e -1 / 3 ), and observe that due to the continuity properties 



limP 



n 



C(s) 



r\{C>n} 



(3.14) 



= 0. 



Finally, letting e tend to zero in (3.13) gives that the probability (3.6) is zero. 

This shows that local maxima are separated with probablity one. In particular, constant 
pieces and no accumulation points of maxima. This concludes its proof. 

Corollary 3.3: Suppose %(i) satisfies the conditions in Proposition 3.1. Then x(t) has 
a.s. only one global maximum in any interval [s, s + 1], t < T. 

To see that Proposition 3.1 is relevant for our problem, we will next show that the process 
3jv(i?) converges to a process of the type covered by this proposition. In fact we have 

Proposition 3.4: Let g : R — > IR + ; g G C°° be an aperiodic even function. Suppose also 
that XiW> $ ^ [0> 27r] is the stochastic process given by 



XiW = g(r(iCOs($- fa)) , 



(3.15) 
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where r is a positive constant, {CijieN; {0i}ieN are two mutually independent families ofi.i.d. 

2 

random variables, distributed as cxe~ x (Q), and uniformly (4>i). Then the process t]n given 
by 

1 N 

rj N (-d) = " E XiW} (3-16) 

v N i= i 

converges in distribution to a strictly stationary Gaussian process r/(i?) having a.s. continu- 
ously differ entiable sample paths. Furthermore, n{jD) has a.s. only one global maximum on 
any interval [s, s + 1], t < it 

Remark: We will use this proposition of course with g(-) = lncosh(/3-). Then the proposition 
implies that the process 3jv(i?) — E<7jv(i?) converges to a Gaussian process with the above 
properties. 

Proof: As are i.i.d. stationary processes on the circle which are infinitely differentiable, 
the convergence of the process to a stationary Gaussian process on the circle is a simple 
application of the central limit theorem in Banach spaces (see e.g. [LT]). A computation 
shows that the covariance of the limiting process is given by 

f(s,t) = E [(xi(s) - Exi(s)) (Xi(t) - Exi(t))] 

2 ( 3 - 17 ) 

= E [g (rd cos(v?i)) g (r& cos(t - s - y?i))] - (E [g «i cos(^i))]) 
We see that this function is even, and is in C°° as a function of r = t — s. Moreover, it 
is easily checked that the covariance function f(r) is strictly smaller than /(0), whenever 
r kir. Proposition 3.1 and Corollary 3.3 then imply the assertions about continuity and 
non-existence of more than one global maximum. This concludes the proof of Proposition 
3.4. 

We now check some intuitive properties of the position of the minimum of the Gaussian 
process from Proposition 3.1 (for those uo such that the minimum exists and is unique). 

Proposition 3.5: Suppose that the conditions of Proposition 3.1 are satisfied. Define 
(jy,jr',P') to be the restriction of (fi, J 7 , P) to all uj such that the conclusions of Proposition 
3.1 are true. Then the position of the minimum 

0*M = arg min XM(0) (3.18) 

i?e[0,7r) 

of the sample path xM is a random variable with uniform distribution on [0, n). 

Proof: To prove that is a random variable, it is enough to show that for all intervals 

U = (a,b) C [0, 7r), the set ■&*~ 1 (U) is in J 7 '. We note that by the continuity of x on [0) 71 ") 
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for all to G 

-d*~ l {U) = {to G SI : x[w](-) assumes its minimum in 

(3.19) 

= {wGO':]teWnQ such that Vs G W n Q, < xOO}- 
The second line can be written as 

|J fl {lo G W : X (t) < X (s)}, (3.20) 

which clearly is in T' . 

Equation (3.20), together with the strict stationarity (since it is a real stationary process) 
of the process %■> implies the uniformity of the distribution. This proves Proposition 3.5. 

Finally, to get some information about the convergence of functions of the position of the 
minimum, we use the following two results. 

Lemma 3.6: Let V([0, n)) be the space of T -periodic, continuous functions, having only 
one global minimum, together with the supremum norm. Then the position $* of the global 
minimum is a continuous function from V([0,ir)) to [0, ir). 

Proof: Suppose that there exists a sequence of functions {f n } converging to / G V([0,n)), 
such that the sequence of the global minima i?* does not converge to •&* , the global minimum 
of /. Then there exists an e > and a subsequence {f nk }, such that for all n^, —■&*] > e. 

Now, since #* is the unique global minimum of /, 35 £ > such that 

f(K k )>f(#*)+5 £ . (3.21) 

Similarly, since is the unique minimum of f nk , 35' e n > such that 

fn k (#*)>fn k W k )+8' e , nk . (3.22) 

Furthermore, since f Uk converges in the supremum norm, V<5 > 0, 3K$ G N such that 

We [0,ir),Vk > K S , \fn h (&) - f(&)\ < S. (3.23) 

For any k > K$ one can therefore write 

fn k m - /on = - umj + /„ fc «j - /«j + /«j - /on 

>S' £ynk -S + 6 £ (3.24) 
>6 £ - 6. 
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Now choose S = ^S £ . Then for all k > Ki s , 

f nk (#*)-f(#*)>l$e>6, (3.25) 

which contradicts the assumption of uniform convergence. 

The following result is crucial to link the weak convergence of the process <7jv($) to the 
weak convergence of the measures Qn,p- 

Proposition 3.7: Define the random sets 

L N [u] = {$ G [0,vr) : r) N [cj]($) ~ min^p') < e N } (3.26) 
with Sn some sequence converging to zero. Then 

L N Z r (3.27) 



Proof: Using the method of a single probability space (see [Shi], Chapter 3, Section 8, 
Theorem 1) one can construct a probability space (fi*, J 7 * ,F*) and random processes rf N , n* , 
such that 

r]* N ^r]*, F*-a.s., (3.28) 

and 

V* = V, V*n = VN- (3.29) 
Now introduce the random level sets 

L* N [u*] = € [0,tt) : V * n [^W)-^V*n[^W) < ?n}, 

Then Ln and L* N have the same distribution. But since rf N [uj\ converges almost surely to 
rf[uj\ G V([0,ir)), one sees that due to Lemma 3.6 £jvM converges P*-a.s. to the position 
of the unique absolute minimum of 77* But this minimum has the same distribution as 
that of 77, which is the uniform distribution by Proposition 3.5. Therefore, Ln converges in 
distribution to a uniformly distributed point on [0, tt). 

We have finally all tools available to prove Theorem 3. 



Proof of Theorem 3: 

F : M(R 2 ) -► R 



We have to check convergence on the following type of functions 
F(p) = F(p(f 1 ),... ,fi(f k )), (3.30) 
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where F is a polynomial function, and fi, ■ ■ ■ , fk are bounded continuous functions from 
M 2 — > R. Convergence in law then means that 



lim E 

JV|oo 



f{QnAA) 



\ r 1 

/ F( — 5l m * cos i?,m* 



1 



sini?) "I - 2 ^(m* cos i?+7r,m* sin $+7r) 



7r) )^. (3.31) 



The left-hand side of (3.31) is explicitly written as 



lim E 

N]~oo 



F{Q N ^[uS\{f 1 ),... ,Q N ,Mih)) 



(3.32) 



We now treat the individual arguments of F in (3.32). Let An [uj] (the level sets in the previous 
lemmata) be decomposed into its 21' connected components Anj n [uj]. As a consequence of 
Lemma 3.7, there exists N[u] which is finite a.s. such that for all N > N(uj), I = 1, and 
the two corresponding connected components are symmetric with respect to the origin. Now 
choose arbitrary points Xnj n [oj] G Anj n [oj]. Then we can decompose 



QjV,/?M(/i) = ^2fi(x NjN )Q Nt f3[oj}(l AN:JN ) + ^2 QNA 1 A N , jN (fi(. x N,j N ) ~ fi)) 



JN 



JN 



(3.33) 



Expanding F using the decomposition (3.33), we get a sum consisting of two different types 
of terms: (i), summands that are products of the first sum on the right-hand side of (3.33) 
only, and (ii) , summands where at least one of the second and third term from the right-hand 
side of (3.33) enter. Proposition 2.3 and Proposition 3.7, and the continuity and boundedness 
of the /j's imply that the terms of type (ii) vanish P-a.s., as N f oo. In the limit, the only 
terms left are of type (i), which together sum up to 



F \ Y,fi( x N, JN )QNA^A N , m ),... ^^nj^QnA^a^J) (3.34) 



JN 



JN 



All arguments of F in (3.34) converge in distribution to 
1 1 

-fi({m* cos m* sini?)) + -fi((m* cos i? + tt, m* sini9 + tt)), Vz = 1, . . . , k (3.35) 

where i? is a uniformly distributed r.v. on [0,tt), by Proposition 3.7. But convergence in 
distribution means by definition that 



limE 

JVToo 
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F(-fi((m* cos i?, m* sini?)) + -fi{(m* cos ~& + tt, m* sint? + tt)) dd, 



TT 



(3.36) 



Section 4 



27 



which in turn is by definition equal to 




m* cos i?,m* sin 




(3.37) 



This proves the convergence in law (1.13) in Theorem 3. To obtain the identification of 
the metastate, just note that the process 77jv(i?)[w] actually converges to the same Gaussian 
process under any of the conditional laws P[-|.F n ], where J- n is the sigma-algebra generated 
by the random variables < n. 00 



4. Volume dependence, empirical metastates, superstates 

We conclude this paper with the discussion of some more sophisticated concepts that have 
been proposed by Newman and Stein [NS2] and Bovier and Gayrard [BG3] and that should 
capture in more detail the actual asymptotic volume dependence of the Gibbs measures. In 
fact, the first question one may ask is whether for a fixed realization as the volume grows the 
finite volume Gibbs states really explore all the possibilities in the support of the metastate. 
One way of stating that this is the case is the following 

Theorem 4.1: There exist (deterministic) sequences N/. j oo such that the empirical 
metastate 



converges almost surely to the law of Qoo,/3- 

Proof: We have seen that the measure QN k ,p is sharply concentrated on the circle of radius 
r* and at the angle where the process g7v fc ($) (defined in (2.3) takes its absolute minimum. 
The idea is to choose Nf, in such a way that these angles will be virtually independent for 
different k. Now note that we can write 




(4.1) 



(4.2) 



where 




i=JV fc -i+l 



(4.3) 
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are independent for different k by construction and 



1 Wfc-i 

Rk(#) = -rrE ln cosh(/3(r*C, cos(tf - (4.4) 



i=i 



Now by standard estimates identical to those presented in Section 3, one shows easily that 
there is a constant C < oo such that 



Nk-i 



<Cexp(-x 2 /C). (4.5) 



Thus we can always choose N k growing sufficiently rapidly (e.g. Nf. = k\) such that R k is 
totally negligible compared to g k for large k, and the position of the absolute minimum of 
<7jv fc (#) is asymptotically equal to that of <?&($). This allows us to approximate for large k the 
random measures Sq n by independent measures and from this the asserted result follows 
from the law of large numbers. 

Remark: Theorem 4.1 says that that the empirical metastate constructed with sparse sub- 
sequences converges to the Aizenman-Wehr metastate, a.s.. This is a special example of 
a general theorem due to Newman and Stein [NS2] (where however they require possibly 
subsequences in the definition (4.1)). 

Rather than considering the empirical metastate with sparse subsequences one may be 
interested in the volume dependence as the volume grows at its natural pace. To capture 
this, the idea put forward in [BG3] is to construct a measure valued stochastic process 

f^ = Km Hp,[tN], (4-6) 

iv Joo 

with t G (0, 1] and to consider either the (conditional) probability distribution of this pro- 
cess (the "superstate" [BG3]) or the (conditional) empirical distribution of the process (the 
"empirical metastate" [NS2]). Let us see what this entails in our context. The reader who 
has been following the exposition of the last two chapters will easily be convinced that this 
problem amounts to study the quantity 

tf(t) = arg mm (xt(e)), (4.7) 

0€[O,7r) 

where Xt(Q) is the distributional limit of the process 



(4.8) 
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where Qn{Q) is defined in (2.3). By completely standard arguments one shows that the 
following invar iance principle holds: 

Lemma 4.2: The process Xtv(^) converges in distribution, as N f oo to the Gaussian 
process Xt{$), t G (0,1], i? € [0, n) with mean zero and covariance 

C(0,^ty) = ^/(<?,0'), (4.9) 



'tt' 
where 

/(#,#') =E[lncosh(/3rCicos(^))lncosh(/3rCicos(v9- (■&-#')))] ■ (4.10) 



XtiP) is a rather curious Gaussian process: as a function of t, to fixed $ it is (normalized) 
Brownian motion, while for fixed t as a function of i? it is the (7°° process discussed in the 
previous section. The question is then what can be said about the process i? t , defined by 
(4.7)? 

Some facts follow easily. For instance, the process is almost surely single valued for all 
t G (0, 1] except possibly on some Cantor set of zero Lebesgue measure. On the other hand, 
it seems natural that such an exceptional set will exist and that a typical realization will have 
continuous pieces and "jumps". Also, for t going to zero, the process "circles" around rapidly 
since \t and \s become uncorrelated as s { 0. But otherwise we do not see any immediate 
more specific characterization of the process or its path-properties. 
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